Specific heat of the skyrmion lattice phase and field-induced tricritical point in MnSi 
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We report high-precision measurements of the temperature and magnetic field dependence of 
the specific heat, C(T,H), across the magnetic phase diagram of MnSi. Clear anomalies establish 
the skyrmion lattice unambiguously as a thermodynamic phase. The evolution of the specific heat 
anomalies, the field dependence of the entropy released at the phase transitions, and the temperature 
versus field dependence of crossover lines provide striking evidence of a tricritical point at {ioHtcp = 
340 mT and Ttcp = 28.5 K. The existence of this tricritical point represents strong support of a 
helimagnetic Brazovskii transition, i.e., a fluctuation-induced first order transition at H = 0. 

PACS numbers: 75.40.-s, 74.40.-n, 75.10.Lp, 75.25.-j 



In the B20 transition metal compound MnSi the para- 
magnetic to helimagnetic transition and the discovery of 
a skyrmion lattice have attracted great scientific interest. 
A comprehensive neutron scattering study recently pro- 
vided strong evidence that the helimagnetic transition 
is driven first order [IH2] by the interactions between 
isotropic critical fluctuations [4], as predicted theoreti- 
cally by Brazovskii [5]. However, the weak first order 
character of the helimagnetic transition has also been 
attributed to the formation of a skyrmion liquid that 
preempts generically the transition to long range mag- 
netic order in chiral hclimagnets [6 8 . A stringent test of 
the Brazovskii-scenario concerns the influence of a mag- 
netic field which is predicted to quench the interactions 
between the critical fluctuations, thereby changing the 
transition from first to second order at a tricritical point. 

The skyrmion lattice consists of a regular arrangement 
of topologically non-trivial whirls in the magnetic tex- 
ture. Its formation has been explained in terms of a 
single thermodynamic phase stabilized by thermal Gaus- 
sian fluctuations [SHU]- However, a competing theoreti- 
cal study predicted a complex generic phase diagram of 
B20 hclimagnets supporting various topological textures 
and meso-phases, driven by a combination of magnetic 
anisotropies and a softened modulus of the magnetiza- 
tion |12j . Putative experimental evidence of this complex 
phase diagram was inferred for FeGe, another B20 com- 
pound, from measurements of the ac susceptibility [T3] 
and essentially featureless specific heat data [Hj. The 
experimental evidence for the complex phase diagram 
was questioned by comprehensive simultaneous measure- 
ments of the ac susceptibility and magnetization in MnSi, 
clearly suggesting the formation of a single skyrmion lat- 
tice phase |15j . The perhaps most important outstand- 
ing experimental challenge concerns thereby evidence for 
clear specific heat anomalies as proof of the skyrmion 
lattice as a thermodynamic phase. 

The precision and density of the specific heat data re- 
quired to resolve these questions does not permit to use 



a conventional heat pulse technique for at least two rea- 
sons. First, typical pulse sizes are of the order of percent 
of the sample temperature. However, to study the ques- 
tions of interest here pulses of at most a few tenths of a 
percent are acceptable to prevent averaging out details. 
Second, the time required for a single data point of at 
least a few minutes renders measurements at very finely 
spaced temperatures and fields unfeasible. 

To obtain high-resolution data for a large number of 
magnetic field values across the entire phase diagram 
we used a quasi-adiabatic large pulse technique on a 
Quantum Design physical properties measurement sys- 
tem (PPMS). This method represents essentially a type 
of adiabatic scanning calorimetry |16j . Typical heat 
pulses had a size of 30% of the sample temperature. 
The temperature was logged as a function of time both 
while heating at a constant power and while cooling after 
switching off the heater. The specific heat reported here 
was inferred from the cooling curves where we applied a 
moving average over 3% of the data points in a relaxation 
curve. The first 1% of every pulse was excluded from 
the analysis. The results obtained are perfectly consis- 
tent with data reported in Ref. [17] where a conventional 
small heat pulse method with pulses of 1% to 2% was 
used. Control measurements with tiny heat pulses of 
0.2 % confirmed that the internal time constants of the 
sample were much smaller than the changes of sample 
temperature. Using this method data at 86 different field 
values were recorded starting from two or three different 
temperatures at each field to confirm the results. It is in- 
teresting to note that temperature scans with a compara- 
ble resolution using a conventional heat pulse technique 
would require nearly one day per magnetic field. 

For our study a large single crystal was grown by opti- 
cal float-zoning under ultra-high vacuum compatible con- 
ditions (for details see Refs. [HI EI]). The single crystal 
was oriented by means of Laue x-ray diffraction and sam- 
ples were cut with a wire saw. The residual resistivity ra- 
tio of samples from the same ingot was 70, characteristic 
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FIG. 1: (Color online) Specific heat of MnSi as a function of temperature for magnetic fields up to 0.7 T along the (110) axis. 
At H = the narrow peak at the helimagnetic transition, marked T c , resides on top of a broad shoulder. In small fields the 
shoulder is characterized by a Vollhardt invariance at T2 . In the field range of the skyrmion lattice the peak at T c splits into 
two peaks where we mark the lower temperature peak as Tax- ( a ) Data as measured, (b) C /T versus T 2 to illustrate the 
lattice subtraction. Panels (c) through (u): specific heat over temperature after subtraction of lattice contributions. Solid lines 
represent polynomial estimates of the background at temperatures far away from the anomalies of the transition. 



of the sample quality studied in most previous investiga- 
tions. The specific heat was measured on a sample with 
the shape of a parallelepiped of roughly 4 x 1.5 x 1 mm 3 , 
for magnetic field parallel to the short side correspond- 
ing to a (110) axis (the soft and hard magnetic axes of 
MnSi are (111) and (100), respectively). For corrections 
of demagnetizing fields (cf. Fig. [2]) the sample was ap- 
proximated as a cuboid J19I . 

To support the specific heat data we also measured 
the ac susceptibility of the same sample for the same 
orientation with a PPMS at an excitation frequency of 
911 Hz and an excitation amplitude of 1 mT. This way we 
avoided systematic differences when comparing the mag- 
netic phase diagrams inferred from the specific heat and 
ac susceptibility, respectively, e.g., due to small differ- 
ences of demagnetizing fields. The ac susceptibility data 
was in excellent agreement with previous work [15 . How- 
ever, to account for a small systematic offset between the 
sample thermometers the temperature scale of the sus- 
ceptibility was shifted by —0.1K. 

Shown in Fig. [11(a) are typical specific heat data of 
MnSi in the vicinity of the helimagnetic phase transition 
as a function of temperature for a wide range of mag- 
netic fields. Our data are in excellent agreement with 
previous reports for the field values available [17J. For 
zero magnetic field a sharp peak at the onset of heli- 
magnetic order, marked as T c , resides on top of a broad 
shoulder. For small fields H, the latter is characterized 
by a Vollhardt invariance j^Oj at a temperature T 2 where 
dC/dH\ H =o = [H US]. With increasing field a sec- 



ond peak emerges in the range of the skyrmion lattice, 
marked as Tai in Fig. [I] This second peak had not been 
reported before. For even larger fields the peak at T c 
shifts to lower temperatures and changes its character 
from first to second order. Finally, we did not observe ev- 
idence suggesting the formation of meso-phases or other 
complexities anywhere. 

The phonon contribution to the specific heat may be 
inferred from the cubic temperature dependence of C(T) 
in the range T 2 < T < 50 K (Fig. 0(b)). The corre- 
sponding Debye temperature, 0d = 515 K, is in excel- 
lent agreement with previous work |17j [21] . Subtracting 
this lattice contribution The detailed field dependence of 
the electronic contribution to the specific heat divided by 
temperature, C e \/T, close to the helimagnetic transition 
is illustrated in the remaining panels of Fig. [I] 

In small applied fields the shape of the peak at T c 
is well described as a symmetric Gaussian with typi- 
cal values of the full-width at half-maximum of AT C ~ 
(0.15 ± 0.02) K. This is characteristic of the latent heat 
of a slightly broadened first order transition where we 
attribute the broadening to imperfections of the sample. 
The peak resides as an additional feature on a broad 
shoulder (Figs, [jjc) through (h)). Up to 140 mT the 
peak at T c is essentially unchanged, while its position 
shifts slightly to lower temperatures. 

Around 145 mT a second peak appears at a tempera- 
ture Tai- This peak clearly shifts to lower temperatures 
with increasing fields and decreases in size until it van- 
ishes above ~ 230 mT (Figs. fl](i) through (n)). In the 
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FIG. 2: (Color online) Magnetic phase diagram of MnSi for 
field along (110) as derived from the specific heat (diamonds), 
as well as the real part (squares) and imaginary part (tri- 
angles) of the ac susceptibility, cf . Ref . |15| . We distinguish 
the following regimes: helical order, conical order, skyrmion 
lattice, fluctuation-disordered (FD), paramagnetic (PM), and 
field-polarized (FP). Phase transitions and crossover lines are 
marked as solid and dashed lines, respectively. In the shaded 
coexistence regime between the conical and skyrmion lattice 
phase a finite Imx(w) is observed, characteristic of hysteresis 
associated with a first order transition. Panel (b) shows an 
enlarged section of panel (a) . 



field range for which two peaks are observed we denote 
the maximum of the peak at the higher temperature as 
Ta2- The shape and size of the peak at Ta2 as well as 
the values of Ta.2 are essentially unchanged in this field 
range. The broadening around 145 mT is due to the ini- 
tial overlap of the two peaks which are well described at 
all fields by Gaussians of widths similar to AT C . 

Finally, for field values exceeding the regime of the 
skyrmion lattice, the shape of the specific heat anomaly 
changes drastically as in Figs.[l](o) through (u) (here the 
maximum is again denoted as T c ). For a field around 
400 mT the anomaly is clearly asymmetric characteris- 
tic of a mean-field lambda anomaly. At the same time 
the broad shoulder on which the peak resides moves to 
lower temperatures. For even larger fields the lambda 
anomaly displays some rounding as shown in Fig. fl](u) 
which may be attributed to the field dependence of T c , 
i.e., temperature scans cut the phase boundary under a 
shallow angle. 

Shown in Fig. pKa) is the phase diagram as inferred 
from our specific heat and susceptibility data. For small 
magnetic fields, the specific heat displays a crossover with 
decreasing temperature from an uncorrelated paramag- 
netic (PM) to a fluctuation disordered (FD) regime. The 
FD regime is characterized by strongly interacting chi- 
ral fluctuations that eventually induce a first-order Bra- 
zovskii transition to the magnetically ordered state [I]. 
This interaction suppresses the correlation length and 
thus induces a point of inflection in the susceptibility, 
d 2 x/dT 2 = 0. Its position defines the crossover line be- 
tween the PM and FD regimes which meets the Voll- 
hardt invariance (VI) for H — > 0. At high magnetic 
fields, where the interaction between the chiral param- 



agnons is quenched, the transition to the conical phase 
is of second order and in the XY universality class. Con- 
sequently, there must exist a fluctuation-induced tricrit- 
ical point (TCP) [22] at intermediate fields where the 
transition changes from first to second order. From gen- 
eral scaling arguments, one expects the crossover between 
the PM and FD regime to merge with the phase bound- 
ary at this singular TCP. Note that also the crossover 
line between the PM and the field-polarized (FP) regime 
observed in the susceptibility emanates from the TCP. 
Taken together, this allows to estimate Ttcp ~ 28.5 K 
and fioHTCP ~ 390 mT. The ordinate on the right hand 
side reflects thereby an estimate of the internal field, 
H int = H - NM « H{\ - N X con), where Xco„ is the 
constant susceptibility within the conical phase. At the 
TCP we find Mo#tCP ~ 340 mT - Around the TCP the 
shape of the specific heat anomaly changes from being 
a symmetric Gaussian to a mean-field like asymmetric 
lambda anomaly. A detailed analysis of the correspond- 
ing entropies to be presented below confirms the position 
of the TCP (see also supplementary information at |23j). 

Decreasing the temperature further within the FD 
regime a first-order transition either to helimagnetic, con- 
ical, or the skyrmion lattice phase is observed depending 
on the applied field. In the field range of the skyrmion 
lattice, two first-order anomalies at Tai and Ta2 with 
Tai < Ta2 are found as a function of temperature. This 
establishes the skyrmion lattice unambiguously as a ther- 
modynamic phase. Interestingly, the dissipative part of 
the ac susceptibility, Imx(w), is quite different at these 
two transitions 15j. A comparatively broad coexistence 
regime with a substantial Imx(w) encompasses the Tai 
line (shaded regime between the conical and skyrmion 
lattice phase in Fig. [21(b)), whereas at the fluctuation- 
induced first-order transition at Ta2 no significant con- 
tribution to Imx(w) is observed. The finite dissipation 
derives here from the nucleation process of topologically 
non-trivial skyrmions within the conical phase [24 , which 
eventually triggers the first order transition at Tai- 

Further, considering the relationship of Clausius- 
Clapeyron, dT/d(fi H) = — AS/ AM , the negative slope 
of Tai(H) at the lower boundary of the skyrmion lattice 
phase implies an increase of entropy as AM > [15] . In 
contrast, there is no difference of entropy when exiting 
the skyrmion lattice phase at larger fields, as the slope is 
vanishingly small. Therefore, the skyrmion lattice phase 
yields a larger entropy than the conical phase. This con- 
jecture is supported by the detailed phase boundary be- 
tween the ordered phases and the fluctuation disordered 
regime shown in Fig. [21(b), where Ta2 is shifted towards 
higher temperatures as expected of a thermodynamically 
ordered state with a larger entropy than the (competing) 
conical phase. 

After subtracting a background to C c \/T determined 
by a polynomial fit (red lines in Fig. [2]) the shaded areas 
were integrated to determine the entropy released at the 




• AS ,o 
■ AS a: 

* AS, 



Na 



.TCP 



% 



I 




0.1 0.2 0.3 0.4 0.5 0.15 0.20 0.25 
^ H (T) n Q H (T) 

FIG. 3: (Color online) Entropy, AStot, inferred from the spe- 
cific heat anomalies at T c , Tai, and Tai- Error bars reflect 
conservative estimates of uncertainties of the background sub- 
traction, (a) For small fields AStot essentially corresponds to 
the latent heat decreasing linearly. Near the TCP critical fluc- 
tuations contribute to AStot resulting in a local maximum. 
The dashed line indicates an estimate of the latent heat that 
vanishes at /xo//tcp ~ 390 mT. (b) Expanded view of the 
regime of the skyrmion lattice. 



phase transitions, AS. The resulting field dependence of 
AS is summarized in Fig. [3](a) where Fig. [3](b) shows 
the field range of the skyrmion lattice in further detail. 
The error bars are estimated by a systematic variation of 
the background contribution; the conclusions discussed 
in the following are not sensitive to the precise choice 
of background (see also the supplementary information 
at [53]). For fields exceeding ^300 mT the anomaly be- 
comes broader and the absolute size of AS depends more 
sensitively on the choice of the background but the gen- 
eral evolution of AS remains unaffected. We denote the 
entropy released near T c by AS to t- In the field range 
of the skyrmion lattice the entropies of the anomalies at 
Tai and Ta2 are denoted as ASai and ASa2, respec- 
tively, where AStot is determined as the sum of ASai 
and ASa2- Note that AS is basically the latent heat for 
the transitions at small fields, whereas AStot corresponds 
to the entropy released near the second-order transition 
for H > Htcp ■ 

The latent heat determined this way is in quantitative 
agreement with Clausius-Clapeyron, suggesting that our 
data capture C(T) without loss of information. On the 
one hand, we find for the lower boundary of the skyrmion 
lattice phase dT / d{n$H) = —(12 ± 1) K/T, correspond- 
ing to dT/d(n H int ) = -(14 ± 1) K/T after correcting 
demagnetizing fields. On the other hand we observe 
around the low field and high temperature boundary of 
the skyrmion lattice phase a change of entropy ASai = 
(1.5 ± 0.2) mJmol _1 K _1 . From magnetization measure- 
ments we estimate for the change at the first-order tran- 
sition AM = (3.4±l)-10- 3 Ai B f-u.^ 1 = (1.3±0.4)kA/m. 
This implies a value of -AM/AS A i = -(13±4) K/T. In 
comparison, at the upper field boundary AM is similar 
to the lower boundary while ASai - ► 0, consistent with 
the absence of a noticeable slope. 



The field dependence of the entropy released at the 
phase boundaries shown in Fig. [3] yields several im- 
portant pieces of information. For small fields, AStot 
can be identified with the latent heat of the first-order 
transitions and it decreases monotonically from its zero- 
field value, 6.8mJmol _1 K _1 , across all phases up to 
~ 300 mT. Second, the unchanged monotonic decrease 
across the skyrmion lattice phase implies that its entropy 
is larger than that of the conical phase as conjectured 
above. Moreover, the field dependence originates from 
the transition between the skyrmion lattice and the coni- 
cal state at Tai, while ASa2 is essentially unchanged. As 
ASai vanishes with increasing field the skyrmion lattice 
phase becomes thermodynamically unfavourable above 
Ha2- Third, near the TCP critical fluctuations start to 
contribute to AStot resulting first in an increase as a 
function of H. In a mean-field approximation, the latent 
heat is expected to vanish as ASi a t cn thoat ~ VHtcp — H 
for H — » if tcp (dashed line in Fig.[3ja)) which is consis- 
tent with the value of -Htcp estimated with the help of 
the crossover lines. For larger fields H > Htcp > the fluc- 
tuations become less singular and AStot again decreases 
until it is finally expected to vanish as the conical phase 
is suppressed by magnetic field. 

In conclusion, we find clear evidence that the skyrmion 
lattice in MnSi represents a thermodynamic phase with- 
out signs of further complexities such as mesophases in 
high-precision measurements of the specific heat. More- 
over, we identify the existence of a tricritical point as 
a stringent test of the hclimagnetic Brazovskii-transition 
at H = 0, i.e., a fluctuation-induced first order transition 
in MnSi. 
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Supplementary Material for: Specific heat of the 

skyrmion lattice phase and field-induced tricritical 

point in MnSi 

In this supplement we present a detailed account of the 
method used to analyze our specific heat data. This con- 
cerns in particular the background we subtracted to de- 
termine the entropy released at the phase transitions and 
the change of its temperature dependence with increasing 
field. 

In the main text we present specific heat measurements 
of the magnetic phase diagram of MnSi in the vicinity 
of the paramagnetic to helimagnetic transition for mag- 
netic fields between zero and noH = 700 mT, see Fig. 

1 of the main text. In order to analyze the evolution 
of the anomalies at the phase transitions determined at 
first a magnetic field-dependent background. Our choice 
of background was defined in the following way as il- 
lustrated in Figs. 1 and 2 of this supplement. We first 
excluded data points in a temperature range close to the 
anomalies. The boundaries of this temperature range are 
illustrated by vertical lines in Figs. 1 and 2. The remain- 
ing data were fitted with a 9 th order polynomial which 
defines a background for a specific choice of boundaries. 
We considered two extreme choices for these boundaries 
of excluded data points. On the one hand, the inter- 
val was chosen as narrow as possible ensuring that the 
background still remains smooth, i.e., without generating 
points of inflection within the background in the temper- 
ature range between the two boundaries. In Figs. 1 and 

2 this is referred to as high background (green line). On 
the other hand, the interval was chosen as wide as possi- 
ble such that the background would still follow the data 
points as closely up to the anomaly as possible. In Figs. 1 
and 2 this is referred to as low background (orange line) . 
As a compromise between these two extremes a middle 
background has been defined (red line). 

In Figs. 1 and 2 data are shown at H = and for typ- 
ical field values in the magnetic phase diagram, notably 
Li H = 180 mT, 280 mT and 390 mT, i.e., a field in the 
skyrmion lattice phase, a field just above the skyrmion 
lattice phase and, finally, a field close to the tricritical 
point, respectively. After subtracting the three different 
background levels defined above, the resulting anomalies 
were integrated numerically to determine the change of 
entropy AS at the transitions. The result is shown in 
Fig. 3 of the main text; error bars represent the two 
limits for the choice of background. 

Regardless of the background subtracted from our data 
the specific heat anomalies display a change of shape as 
a function of temperature. This is illustrated in Fig. 2, 
where the anomalies have been fitted by Gaussians (thin 
blue line), 

AC , (T-T c ) 2 

— =y + Ae ~^^- (1) 

Gauss 



The parameters of the fits, A, w and T c , are shown 
in each panel. The offset yo is below O.lmJmol -1 K 2 
and thus negligible. At zero field and low fields, the 
anomalies are symmetric and essentially Gaussian; at 
high fields the shape cannot be fitted by a Gaussian. 



This is shown for the field fioH = 390 mT close to the 
tricritical point which clearly demonstrates an incipient 
asymmetric lambda anomaly characteristic for a second 
order transition. 
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FIG. 1: Typical specific heat anomalies and step-by-step analysis to determine the change of entropy at the anomalies as shown 
in Fig. 3 of the main text. For close-up views see Fig. 2 in this supplement. Columns correspond to the same field value. The 
top row shows the data as measured after subtracting a phonon contribution; colored vertical lines indicate the temperature 
interval for which data was excluded when determining the background (red line). The centre row shows the specific heat 
anomalies after subtracting different choices of background as explained in the text. The bottom row displays the entropy as 
inferred from the specific heat anomaly by numerical integration. Data points in Fig. 3 of the main text were determined from 
the non-zero plateau values; the variation of these values with the choice of background defines the error bars in Fig. 3 of the 
main text. Panels (af) through (cf): behavior at zero field. Panels (a2) through (c2): behavior at fj,oH — f80 mT, i.e., a field 
value in the skyrmion lattice phase. Panels (a3) through (c3): behavior at fj,oH = 280 mT, i.e., a field value just above the 
field range of the skyrmion lattice phase but below the tricritical point. Panels (a4) through (c4): behavior at noH = 390, i.e., 
close to the tricritical point. 
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FIG. 2: Close-up view of the data shown in Fig. 1 of this supplement. The specific heat anomalies at low fields are well 
described by a Gaussian but turn asymmetric at high fields indicating the presence of a tricritical point around /io-ff > 390 mT. 
Columns correspond to the same magnetic field value. The first row shows the data as measured after subtracting a phonon 
contribution. The second, third, and fourth row show the specific heat anomaly after subtracting the high, middle and low 
background, respectively. The thin blue lines represent Gaussian fits to guide the eye. The shape of the anomaly changes 
with increasing field from a symmetric, essentially Gaussian shape to an asymmetric shape, which clearly cannot be fitted by 
a Gaussian. This conclusion is independent of the background subtracted. 



